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Abstract 

Mutual equilibrium in long-range interacting systems which involve nonadditive energy, is ef- 
fectively described in terms of entropy with a nonadditive composition rule. As an example, long 
range Ising model is considered. The generality of the term having product of the system entropies 
is pointed out in this framework. 
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1. Introduction 



3e the correct entropy of the system, 



Additivity of quantities such as entropy and energy, is an essential premise of thermody- 
namics and statistical mechanics of systems with short range interactions [lj. Arguments 
about positivity of specific heat or concavity of the entropy-energy curve are based on this 
additivity postulate. On the other hand, systems where long range interactions are relevant, 
are basically nonadditive. Due to this feature, thermodynamics of such systems displays 
unusual properties, like inequivalence of different ensembles, negative specific heat in micro- 
canonical ensemble and possible temperature discontinuity at first order transitions |2|. For 
nonadditive systems, the non-concave entropy ma> 
with the result, the specific heat may be negative ;c 

Similarly, the notion of macroscopic thermal equilibrium, which introduces the concept 
of temperature through the zeroth law, is based on additivity of entropy and energy. On the 
other hand, for large systems with long range interactions or even finite systems with short 
range interactions, the concept of equal temperatures over subsystems, does not hold at 
thermal equilibrium (maximum of total entropy). Clearly, if additivity postulate is relaxed 
for either entropy or energy, we have breakdown of the zeroth law. 

In this paper, we consider two macroscopic systems in thermal contact with each other, 
whose particles experience long range interactions. The total system is isolated and so its 
energy is fixed. We consider the problem of maximization of the total entropy under the 
constraint of fixed total energy. Two alternate approaches are explained: the first one is 
based on the locality and local extensivity assumptions 6]. Here all the correlations due 
to interactions are completely described by correlations in the energies of the systems. The 
total energy E is nonadditive and in general given by 

E = E 1 + E 2 + G(E 1 ,E 2 ), (1) 

where E\ and E 2 are the energies of system 1 and system 2 respectively, in the absence of 
interactions G(Ei, E 2 ). These are called local energies. Note that these are not the current 
energies of the subsystems which become modified due to the presence of interactions. The 
total entropy of the two systems, written in terms of E\ and E 2 , is additive: 

S(E 1 ,E 2 ) = S 1 (E 1 ) + S 2 (E 2 ). (2) 

Now upon maximization of the total entropy under constraint of Eq. (0), one obtains that 



temperatures of the two systems are not equal and the procedure fixes only the ratio of the 
two temperatures, unless entropy is defined by, say Boltzmann's principle. 

In this paper, we also study the approach, complementary to the above, where all correla- 
tions due to interactions are absorbed in the entropy. The entropy so defined will naturally 
obey a generalized (nonadditive) composition rule. On the other hand, the energies in this 
alternate discription will be taken as additive. Thus our motivation is to see how nonad- 
ditive energies are effectively described in terms of the nonadditive entropy and what that 
entropic rule actually looks like. 

The paper is arranged as follows: In section 2A, we treat the case when Eqs. ([Q and (j2J 
are assumed; section 2B presents the alternate approach whereby energy is assumed to be 
additive. Thus we introduce the properties of generalized entropy and do the maximization 
with an additive energy constraint. In section 3 A, we take a specific long-range interacting 
model, namely a spin-lattice model, and describe the maximization of entropy as in the 
approach of section 2A mentioned above. The final section 4 is devoted to conclusions. 

2A. Generic nonadditivity of energy 

Let us consider the variation of the total energy E — E ( E\ , E2 ) , 



From Eq. ([TJ, we obtain 



and 



' dE\ 1 ( dG\ 



dE\ 1 ( dG\ 



Due to constraint of the fixed total energy, we obtain 

dE = dEi = -dE 2 . (6) 

Now if the total entropy is given by eq. (J2J), we also impose dS = 0, which implies 
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or 



dSi dS 2 /0 s 
^dE, = Ql dE 2 (8) 



Here j^h = ^- is defined to be the inverse local temperature of each system p. Clearly, 
these temperatures are not equal in the presence oflong range interactions. 



2B. Second approach: Nonadditive entropy but additive energy 

Now we formulate an equivalent description in terms of nonadditive entropy. We assume 
that all the previous correlations expressed through energy are absorbed in the generalized 
entropy. Thus the energy in our alternate description is taken as additive. Our goal is to 
seek the composition rule for entropy, so that upon maximization of total entropy under the 
constraint of fixed total energy, we obtain a form similar to Eq. (JHJ) • 

We suppose that the generalized entropy of the composite system made up of two systems 
1 and 2, can be written as 

S(l,2) = f(S u S 2 ), (9) 

where Si(Ei) and S 2 (E 2 ) ar e the system entropies, as functions of respective local energies. 
Note that although we use here the same symbols for entropies and energies as in the previous 
section, these two sets, clearly are different quantities. / is a certain bivariate function, such 
that the function itself and all its derivatives are continuous. 

Thus consider the sitation when the total energy of the system is fixed at E. Let Ei 
and E 2 be the energies of the systems when the total entropy /(Si, S 2 ) is maximized. Thus 
putting df = 0, under the constraint of additive total energy; 

E = E 1 + E 2 , (10) 

we obtain 

f df\ dSi ( df\ dS 2 . , 

With dE = 0, we have dE\ = —dE 2 . Thus Eq. (II lj) can be written as 

df_dSi dldS, 

8S 1 dE 1 dS 2 dE 2 1 ] 

Now when we seek a description equivalent to the first approach, we require that Eq. ()12j) 
be of same form as Eq. (jSJ). Thus the following conditions hold 

dG 



1 + 



dE 2 



(13) 



° J <f>(E h E 2 ) 



dG 



(14) 



dS 2 

where <f> is some unknown different iable function. Here it is apparent that the composition 
rule for the generalized entropy will turn out to be ambiguous due to the arbitrary function 
<p. However, it is possible to make a specific choice, as following. Note from Eq. (0) that 
when G = 0, the total energy is additive over local energies of the systems. In this situation, 
it is expected that = 1, i.e. entropy is additive. Thus in the following, we set the function 
4> as a constant equal to unity. 

Now we evaluate the following cross derivatives 

d 2 f d 2 G dE 2 d 2 G 

0S 2 dS 1 8S 2 dE 2 dS 2 dE 2 2 [ ' 

d 2 f = d 2 G dEi d 2 G 

dS 1 dS 2 dS 1 dE 1 dS 1 dE 1 2 ' [ ' 

Assuming that the function / and its derivatives are continuous, we get the equality of the 
above cross derivatives, yielding 

dEi^G_ _ dE^d^ 

'dS^dE^ 2 ~lS 2 "dE^' ( ' 

We simplify the above condition by invoking separation of variables. For this we assume 
that J^rl, involve variables or parameters of the respective system (i) only. The simplest 
case is when 

be? = *> (18) 

where Ui is constant and may depend upon the pre-specified parameters of the system. Thus 
we have 

dE\ dE 2 

ds;" 1 = IT 2 = x ' (19) 

where \ is a constant. These equations can be integrated to give 

Si = -Ei + ai, (20) 

X 

where (i — 1,2) and Oi is constant of integration. 

Moreover, G of Eq. Q is a known function of E\ and E 2 , so that 

BG 

— = u t E t + ip(Ej) (21) 

where i j and if) is a known function. 

Eq. (|2*Tj) can be expressed in terms of Si(j) using pOj). Then Eqs. (fTSjl and (JHJ) (with 
0=1) can be integrated to obtain the generalised composition rule for the entropy. 



3A. Long range Ising model: an example 



In this section, we take a specific long range interacting model as an example of the 
general treatment presented so far. Consider a lattice of N spins with total energy as given 
by 



E = h x e x + h 2 e 2 



Ji 2 ^2 2 t 



(22) 



2 2 

The total system is assumed to be separated into two regions with a fixed number of sites. 
The spin excess = Nf - Nf, where iVf are the number of up (+) and down (-) spins 
inside each region. Jj and hi are the known coupling constants within each region. Ji 2 
represents the coupling between each region, assumed to be smaller than Jj. Note that the 
local energy of each region is Ei = h^i. Writing E in terms of E\ and E 2 , consider the 
variation of the total energy, dE{Ei, E 2 ) and note that 



and 



'BE 
dE 1 



'bet 

dE 2 . 



E 2 



Hi 



1 ^1 P ^12 T7I _ 

1 — -— 9-Ga — ——hj 2 = ax, 
hi h x h 2 



h 
h 2 



—E 2 



J 
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h\h' t 



Ei = a 2 . 



(23) 



(24) 



Imposing the maximisation of total entropy (additive) under the constraint of a fixed total 
energy E, as discussed in section 2 A, the equilibrium condition, Eq. (JHJ) takes the following 
form 



1 



hi 



—E 2 



J 
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hh 



Ei 



dSi 
dE 1 



Jl E 
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h\h 2 



Eo 



dS 2 
dEo 



(25) 



3B. Long range Ising model: second approach 



Following the general procedure, described in the beginning of section 2B, we now demand 
that applying the second approach to this model we obtain the same equilibrium condition 
as ([25)1 . Thus the following conditions hold 



df_ 

OS! 

df_ 

dS 2 



<P(E U E 2 ) 
(f>(Ei,E 2 ) 



1-4^2 

h 2 2 



12 



hih^ 



hi hh 2 



(26) 
(27) 



Here also, as in section 2B, we restrict to the case when 



1. 



Then equality of the cross derivatives of / yields the following condition 

J\ dEi J 2 dE 2 



(28) 



hi 2 dSi h 2 dS 2 

Thus the constant v± of Eq. is expressed in terms of the parameters of the individual 
systems. Now, as each side pertains to an individual sytem, so we equate each side to a 
constant x- Finally, we obtain 

st = W- i = 1 ' 2 ' (29) 

which effectively implies a linear relation between the entropy of the system and its energy 

Si = -^Ei + <n. (30) 

Xhi 

Here cr, is constant of integration. It can be determined from the fact that Si — for 



Ej = —h;Nj. Thus O; 



Using the above relation, we can write Eq. (J2HJ) (with = 1) as follows 



df 



1 - x(s 2 - (r 2 ) - ~ a i) 

J\tl 2 



(31) 



dSi 

and a similar relation corresponding to Eq. f)27j) . These expressions can be easily integrated 
to obtain the form of total entropy in terms of the system entropies 

f(Si, S 2 ) = Cl S x + c 2 S 2 - ^f^S, 2 - ^f^S 2 2 - 2XSA. (32) 

This is a form of nonadditive entropy rule that effectively describes the mutual equilibrium 
between two long range Ising model systems. The constants q are explicitly given by 

Ci = 1 + —— + —r — , 33 

n 2 n 2 

C2 = 1 + — + — (34) 

4. Conclusions 

We have mapped the mutual thermal equilibrium of a bipartite long range interacting 
spin model, in terms of additive total entropy and nonadditive total energy, to a description 
in terms of nonadditive total entropy and additive total energy. The composition rule for the 
nonadditive entropy has been derived. Additionally, we make the following observations: for 



J12 = 0, the total entropy is in the form / = ciSi + C2S2 — 2\SiS2- The additional product 
term is due to the fact that there are still interactions within each system, represented by 
the coupling constants J\ and J 2 . 

We comment on the generality of this product term, which means it always appears in 
the entropy composition rule, when in Eqs. (J 17)) we assume that 

= ^), (35) 

where g is some integrable function of its argument. 

To show this, note that due to separation of variables, we can write 

f^) = f^) = «. < 36 > 

where £ is a constant of separation. These equations can be integrated to give 

£Si = J g(Ei) dEi + const. (37) 

Integrating Eq. (J3H|) and using Eq. IpTTj) in the result so obtained, we can write 

dG 

— = £Si - const + V>(#;) (38) 

where i ^ j and -0 is a known function. Again, using the above Eq. in Eqs. (|13|) and (|14j) . 
we see that the obtained / function after integration, necessarily contains the term ^SiS^. 
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